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Abstract 

We discuss the energy scale profile of the bi-maximal mixing which is given at the GUT 
energy scale in the minimal SUSY model, associated with an assumption that Yj^Y^, is 
diagonal, where Y^y is the neutrino- Yukawa coupling matrix. In this model, the Dirac 
mass matrix which appears in the seesaw neutrino mass matrix is determined by three 
neutrino masses, two relative Majorana phases and three heavy Majorana masses. All 
CP phases are related by two Majorana phases. We show that the requirement that the 
solar mixing angle moves from the maximal mixing at GUT to the observed one as the 
energy scale decreases by the renormalization effect. We discuss the leptogenesis, and 
the lepton flavor violation process by assuming the universal soft breaking terms. 
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1 Introduction 



The bi-maximal mixing scheme[l] may be most attractive one. It has a simple and 
beautiful structure and there are various models which give the bi-maximal mixing at 
the GUT scale[2] . In addition to this, we feel the property Vis = is interesting. If the 
bi-maximal mixing is realized at the GUT scale Mx, the Dirac CP phase 6 as well as |Vi3| 
which are absent at Mx are induced following to the renormalization group equation at 
the low energy. This may give us a chance to predict these quantities. Another interesting 
point is that we may able to solve the discrepancy between the maximal solar mixing 
angle at the GUT scale, tan^6'0 = 1 and the experimental data[3,4] at the low energy 

tan^ - 0.40 . (1) 

Let us consider the renormalization group equation due to the neutrino- Yukawa and 
the T- Yukawa couphngs. We have shown[5,6] that the effect due to the r- Yukawa cou- 
pling rotates the solar angle toward the dark side. Therefore, the large neutrino- Yukawa 
couplings are needed [6, 7] to compensate this and rotates the solar angle toward the nor- 
mal side. In this analysis, Majorana CP phases[8] in the neutrino mixing matrix[9] play 
an important role. 

In this paper, we continue this analysis further by considering the neutrino mass 
matrix derived through the seesaw mechanism in the framework of the MSSM with the 
universal soft supersymmetry breaking terms. Our main motivation is to examine the 
the structure of the Dirac mass matrix and explore the possible relation of CP phases at 
the high and low energies. In general, the Dirac mass matrix introduces new CP phases 
and there is no relation between CP phases in the low energy and the high energy [10]. 

In Sec. 2, we explain the assumptions which we adopt to construct the neutrino mass 
matrix. The renormalization group analysis is briefly explained and a typical form of 
the neutrino- Yukawa couplings is discussed. By assuming this typical from, the Dirac 
mass is determined. In Sec. 3, the various results including the asymmetry parameter 
of the leptogenesis, the lepton flavor violation are discussed. The numerical analysis is 
presented in Sec. 3. In Sec.4, the summary and discussion are given. 
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2 The model 



We assume the MSSM with the neutrino- Yukawa couphng matrix, Y^, and the right- 
handed Majorana neutrino mass matrix, M. The related terms are given by 

1- 



>C^+M = NR(PiY,UL - -{Nr)CMNr + h.c. . (2) 

The soft SUSY breaking terms arc assumed to be universal and the source of the Icpton 
flavor violation (LFV) is only through the neutrino- Yukawa couplings. The left-handed 
neutrino mass matrix, m,^ is derived through the seesaw mechanism as 

m,{Mx) = mlM-^mD , (3) 

where 

run = (4) 

with Vu = V sin j3. Here, the neutrino mass matrix is effectively given at the right-handed 
neutrino mass scale, while Y^, and M in Eq.(3) are defined at the GUT scale, Mx- 
The Dirac mass matrix is generally expressed by 

mD = V^DdVl , (5) 

where Vr and Vl are unitary matrices. Do is a diagonal mass matrix 

Dd = diag(m£,i, mD2, mm) , (6) 

with real, positive eigenvalues, rriDi. In the following, we take the diagonal basis of M, 

M = L>« = diag(Mi,M2,M3) , (7) 

with real positive eigenvalues, Mj. 

In this paper, we consider the hierarchical Dirac mass case, 

rriDs > "^D2 > rnoi ■ (8) 

Also we take 

Mg > M2 > Ml . (9) 



(1) The model 

Our mode consists of the following contents. 

(A.l) The bi-maximal mixing is realized at the GUT scale, Mx- 

(A. 2) The experimental solar mixing angle is achieved by the renormalization group 
effect due to the neutrino- Yukawa and the r- Yukawa couplings. 

(A. 3) Y^Y^ is assumed to be a diagonal matrix. 



Prom (A.l), the neutrino mass matrix, rrii, at Mx is given by 

m,{Mx) = ObD.OI , 

where Ob is the bi-maximal mixing matrix 

-75 M 

^2 2 ^ ' 

Di, is a diagonal matrix with complex eigenvalues. 



(10) 



= diag(mi, m2, ms) = diag(|mi|, |m2|e*"°, |m3|e*'^°) , 
where ao and /3o are Majorana phases[8]. 



(12) 



(2) The motivation of the assumption (A. 3) 

In our papers [5,6] , we showed that the r- Yukawa contribution rotates the solar angle 
toward the dark side. Therefore, (A. 2) requires[6] that the contributions from the 
neutrino- Yukawa should compensate the r- Yukawa effect and rotate the solar angle 
into the normal side. For this, some elements of Y^Y^, must large. On the other 
hand, the LFV processes take place through mixings in the slepton sector [11]. The 
rate for £j —>■ ij+'jis proportional to |(Fjyi,)jj|^, and if we require that its branching 
ratio is less than, say, 10^"*^^, then |(Fj'Y',^).y | < 3 x 10~^(1/ tan/3) {i ^ j) is required. 
Thus, only the elements which can be large are diagonal elements of YjYi,, which 
leads to the the assumption (A. 3). The possible modification of this assumption 
will be mentioned later. 
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(3) The form of mo = {vu/V^)Y^, 

In the following, we sometimes use rriD instead of for the convenience. From the 
assumption (A. 3), is expressed by 

rUD = V^DoPex , (13) 

where Vr is a unitary matrix, Dd is a diagonal matrix defined in Eq.(6) and Pgx is 
the matrix to exchange the eigenvalues oi Djj. 

The matrix Pei is fixed by considering the renormalization group effect to the solar 
neutrino mixing parameters. We take 

/ 1\ 



p 

ex 



1 

V-1 oy 



(14) 



from the reason we explained later. 

(4) The renormalization group effect 

The renormalization group equation is given for Mx > n > by 

^ = {[(>;'^^)^ + (Ylyer\m. + mA{Y}Y^) + (y;i;)]} , (15) 

aside from the terms proportional to the unit matrix. Here, Mr is the right-handed 
neutrino mass scale. For the charged lepton Yukawa coupling matrix, we consider 
only the r- Yukawa couphng, Y^. — diag(0, 0, y^). When ji < Mr, only the r- Yukawa 
couplings contribute, because the heavy neutrinos decouple from the interaction. 

From the assumption (A. 3), we express 

y;y, = diag(i/^, 1/2,1/2). (16) 

The contribution from it is split into the one proportional to the unit matrix, say, 
diag(0, y|,0) and the rest, diag(yi — y|, 0,^/3 — y|). The former contributes the 
overall normalization of neutrino masses so that we discard it. As a result, the 
renormalization equation is expressed in a good approximation as 

m^{mz) = m^{Mx) + K'^m^{Mx) + m^{Mx)K , (17) 



where [6] 

/ee 0\ 
K= 
VO ej 
The Ee and are given exphcitly by 



(18) 



yl - yl 


In 




167r2 


KMr 


yl-yl 


In 


{mx 


167r2 


KMr 



+ J.ln^^^ 



167r2 



(19) 



Here we neglect the threshold effect of Mj, which we discuss later and take Mj = 

Mr. 

Since the renormalization group effect is discussed in detail in Refs.4 and 5, we give 
only the result. The effect to the sizes of neutrino masses, the atmospheric mass 
squared difference and the atmospheric mixing angle are small for \mi\ < O.leV. 
Only the effect appears to the solar mixing angle and the solar mass squared masses, 
which are related by 



tan^ dr: 



1 + (e^ - 2ee) cos2(ao/2)m?/Am| 











(20) 



l-{er- 2ee) cos2(ao/2)mf/Am| ' 

where Oq and Amg are the experimental values of the solar mixing angle and the 
mass squared difference which are defined at the low energy scale, mz- In order to 
obtain tan^ 9q ~ 0.40, it is required 



(2ee - e-,) cos^ {ao / 2) {ml /AttIq) = cos 26*0 



(21) 



This equation gives the constraint on yfml. Therefore, the smaller neutrino masses 
requires the larger Yukawa couplings. Since the Yukawa couplings can not be 
very large, the neutrino masses must be large. If we take mi — O.OSeV, Amg^i = 
6.9 X lO-^eV^ cos 2^0 = 0.43 and | cos(ao/2)| = 0.5, we need ~ 0.5. 

The important point is that the condition 2ee — > is necessary. The r- Yukawa 
couphng gives €e — and > 0, so that this effect rotates the angle into the dark 



side. To compensate the r- Yukawa contribution and rotates it into the normal side, 
we need the large neutrino- Yukawa couplings and they satisfy 

'2yl>yl + yl + yl- (22) 

This is satisfied only when the neutrino- Yukawa coupling matrix has the inverse 
hierarchical structure. That is, we have to assign 

V2 



yi = 

y2 = rnD2 



V sin (3 ' 

V sin (3 ' 



ys = rnm ^ ^ , (23) 
which fixes the form of Pex in Eq.(14). 

(5) The form of Vr 

In the diagonal basis of M = Dr, we have 

ObD.OI = mlD^^mD ■ (24) 

By substituting Eq.(13) and after some computations, we have 

Mr' = {V*^D-^^Vl) = Dl\P,,OB)D,{Pe.OBfD-^'^ . (25) 

We observe that Vr and Mj are determined by complex neutrino masses rrii and real 
Dirac masses moi, ie., 6 real positive masses and two Majorana phases. In other 
words, all CP violation phases in this model are related to two Majorana phases. 

Now we diagonalize M^^ in Eq.(25) under the condition in Eqs.(8) and (9). In 
addition, we assume that 

mi ~ \m2\ ~ \mz\ , (26) 

because rrii should be much larger than \fKm\^^ for the renormalization group to 

be effective. We define 

a = mi m2 2m3 , 



and 



which satisfy 



b — mi + 777,2 — 27773 , 

c = V2(m2 - mi) , (27) 



5i = ^,S.^^, (28) 
mD2 runs 



(5i ~ ^2 < 1 . (29) 



The diagonahzation is exphcitly given in Appendix, so that we give the result. The 
unitary matrix Vr is given by 

a ^ a+b ^ ^ 



R 



Si 1 -^^2 

a+b ^ 



(30) 



where eigenvalues are 



Ml Ami, ' 

e^'f'^ _ 1 4(0 + 6)7713 

M2 47?7|,2 a 

e**^^ 1 16?77i?772 

M3 ~ 4777?,o a + 6 ' 



(31) 



'3 ^"^733 

where phases, 0^ are defined such that Mj are real positive. Here we see in general 
M3 ^ M2 ^ Ml because m/^a ^ r?7D2 ^ ^di, except for the very special case 
777i + 7712 = 0. It may be interesting to observe that 

77iim2m3MiM2M3e-^('^i+'^^+'^3) = ml^m^^^il^ ■ (32) 

As you see in Eqs.(30) and (31), the mixing matrix Vr and also the heavy neutrino 
masses M, arc determined by \mi\, two Majorana phases, i.e., their relative phases 
and moi- Since the Dirac mass matrix is given by 777^ = VrDdPsx, phases in it is 
determined by two Majorana phases. 



3 The numerical analysis 



In the following, we keep only the yi term in the neutrino- Yukawa couplings. Also, we 
take 

mi = \m2\ = m , 

Imsl = m±^^. (33) 

In the following, we take the solar squared mass difference given by KamLAND coUabo- 

ration[12], Amg^j = 6.9 x lO^^(eV)^ and the atmospheric squared mass difference by the 
SuperKamiokande collaboration [13], Aml^^ = 2.5 x 10~^(eV)^. The atmospheric mixing 
angle is the maximal as we chose the bi-maximal mixing scheme. 

(1) The information from the renormalization group analysis 

(1-a) The right-handed neutrino mass scale 

If we consider the heavy neutrino threshold correction, ee in Eq.(19) should be 
modified by 

1 



' 167r2 

where 



(^J>Cy,)n , (34) 



C = diag(ln(Mx/Mi), ln(Mx/M2), ln(Mx/M3)) . (35) 

We obtain 

1 '"D3 ,„ 



For €t-, Us — yi is neghgible because of Eq.(8). By substituting eg and in 

Eq.(19) into Eq.(21), wc have 

mMs , fMx\ 47r2 cos 26'0f 2 Am?, sin^ /3 m^tan^/?, /Mx\ 

In hrr = /o^ + n 111 z— '(37) 



|cos(ao/2)| VMa/ m2cos2(ao/2) 2 \mz 

where we used mf^g ~ mMs/\ cos(q;o/2)|. 



M3[GeV] 



1.2-10 



1-10 



8-10 



6-10 



4-10 
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14 



14 



13 



13 



13 



13 



fm = 0.025eV 
m = 0.05eV 
m = O.leV 



0.5 



0.6 



0.7 



0.8 



0.9 



COS 



Figure 1: The relation between M3 and | cos(ao/2)| for m = 0.025eV, 0.05eV, O.leV. We 
take tan^ Oq = 0.4, tan (3 = 20, and Mx = lO^'^GeV. 
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In Fig.l, we show the relation between M3 and | cos(q;o/2)| for m =0.025eV, 
O.OSeV, O.leV and tan^ 9q = 0.4 with tan/3 = 20, Mx = lO^^GeV. The M3 

behaves almost independently of | cos(q;o/2)| for m = O.OSeV. 

In particular, when |cos(q;o/2)| = 0.50, cos 2^© — 0-43 and m — O.OSeV, we 
get 

M3 ~ 4 X lO^^GeV , (38) 

which implies 



mD3 = \ I ,1 - 60GeV . (39) 

V I cos(q!o/2)| 

If we take tan/3 larger than 20, M3 becomes larger. The other masses Mj(i = 
1,2) are determined once mDi-i olq and (5q are given. 

(1-b) The induced |Vi3| 

The I Vis I and the Dirac CP phase are induced. We show only |Vi3| which is[6] 
W I - 010 ( l^-l ( "^'"^^ \ f 2.S X 10-=^eV^ \ / sin(ao/2) \ 

where the value — 9.2 x 10~^ is the one for tan/3 — 20. Thus, the model 

generally predicts the value of |Vi3| which is consistent with the CHOOZ's 
bound[14] and may be detectable in the near future experiments. 

(1-c) The Dirac phase 5 

The induced Dirac phase is given by[6] 

Y-/3o-| + 6 + 6, (41) 

where ^1 = arg(c — se~*°°/^) and ^2 = arg(c + se'"°/^). Here, c = cos^ and 
s = sin ^ and 



so that 



sin 29 cos(ao/2) = - cos 2^© , (42) 

|cos(q;o/2)| > cos2eo ■ (43) 
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The relation between S + Po and | cos(q;o/2)| is numerically plotted in Fig.l of 
Ref.6. For | cos(ao/2)| ~ cos2^^0 ~ 0.43, 6 + /3o takes values between — 7r/2 
and —3n/2. In the discussion of the leptogenesis, we show /3o ~ is favored to 
reproduce the experimental value of the baryon asymmetry. Then, for Pq — 0, 
our model predicts 

--7r>S> --TT . (44) 

(l-d) The neutrinoless double beta decay 

The effective mass[15] of the neutrinoless double beta decay is given by[6] 

(mix) ~ m| cos(q;o/2)| , (45) 

Since |cos(q;o/2)| > cos26'0, (m^) > mcos2^0 ~ 0.43m, which may be 
within the experimental sensitivity in the near future. Our expectation is 
m ~ 0.05eV, so that it may be around 0.02eV. 

(2) The leptogenesis 

The lepton asymmetry parameter, e is defined by [16] 

^ _ r{Ni M^) - r{Ni ^U) 



= 7^7 \—j:MmDml)yXMf/Ml), (46) 

where v ~ 246GeV and 

(47) 

The approximate form of / is valid for our case because Mj have the hierarchical 
structure as in Eq.(31). 

By the explicit computation, we find 



Dl 



\a\ 


l' + 


\b\' + \ 


|C|2 




\a\ 


2 



b \c 



2 



(mz,ml,)i2 ~ m^im^2 - + , ' ^ J e'^^'-^'^^ , 
(mz,mlj)i3 ~ m^im,,3 (^) e^^^^-'^^)/^ (48) 
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By using Eq.(48), e'"'^^ = (Mi/4m|)i)a*, e'"'^^ = (M2/4m|)2)(4(a + 6)m3/a)* and 
g-i<^3 ^ (M3/4m|53)(16mim2/(a + 6))*, we find 

Now, we find up to the first order of Am|i/m^, 

_ 3mMi cos(ao/2) / Amji \ ap 
' - 167r^2 ) ''""^Y ~ ' ^^^^ 

wliere 



= y 1 + cos2 ^ + 2 cos ^ cos(^ - /3o) , (52) 

and Am|i = Amatm for jmsl > |mi| and — Amatm fo^ I "^3 1 < It may be 

commented tliat in tlie approximation of mi = \m2\ = 1, the contributions from 
Eq.(49) and (50) cancel each other, so that e is suppressed by {Aml^^/m^). 

In Fig. 2, the asymmetry parameter e is plotted as a function of | cos(q;o/2)| and Pq 
with m = O.OSeV, 'm£,i/m£)2 = m/)2/''^D3 = V^j and tan/9 = 20. The larger e is 
obtained for smaller | cos(q!o/2)| > cos26'0 and also Pq. The black lines show the 
numerical computation without any approximation and the gray lines are obtained 
by using our approximate formula in Eq.(50). For | cos(q;o/2)| < 0.6, there are some 
difference between the exact computations and Eq.(51). We can obtain e ~ 10~^ 
for /?o ~ and cos 26^0 < | cos(ao/2)| < 0.8. 

The baryon asymmetry parameter is given by [17] 

rjBo ^ -10~^eKo , (53) 

where for small m, — l/(2/PT9) with K - 170(m/eV). With m = 0.05eV, 
we find ~ 6 x 10~^, so that we obtain 

77BO ~ 6 X 10-1° , (54) 
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Figure 2: The asymmetry parameter e as a function of |cos(a;o/2)| and /Sq. We take 
m = O.OSeV, mDi/mD2 = ?7^d2/'"^D3 = 1/5, and tan/? = 20. The black hues show 
the exact computations and the gray hues show the numerical values obtained by the 
approximate formula in Eq. (51). 
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which agrees with the experimental value [18]. 

Of course, the value of the asymmetry parameter depends on Mi linearly, which 
we derived by assuming mDi/mD2 — '>tid2/'>tT'D3 — 1/5. If the hierarchy of the 
eigenvalues of the Dirac mass matrix is larger than what we used, we find the 
smaller Mi which results in the smaller asymmetry parameter. In other words, we 
can explore how hierarchical the Dirac masses are from the asymmetry parameter 
e. 

(3) The LFV processes 

In this model, the LFV processes take place through the slepton mixing, which is 
absent at the GUT scale. However, the slepton mixing is induced by the renormal- 
ization group effects at the scale Mr where the right-handed Majorana neutrinos 
are decoupled. In the leading log approximation, the off-diagonal terms of the 
slepton mass matrix is given by [18] 

i'^ih - — ^Y^2 iV^Yi^h , (55) 

for i ^ j, where C is given in Eq.(35). The off-diagonal elements contributes to 
lepton flavor violating processes such as — > 67. The decay width of li Ij^ 
process is approximately given by[18] 

^ (6rng + 2Ag)^ tan-^ ^ 

where ms represents typical mass of supersymmetric particles. 
It is convenient to separate m\)£mD into two parts as 

m\)CmD =< > mj^rriD + m}jj{jC— < C >)mD , (57) 

where 

<£> = ln(Mx/M2)diag(l,l,l) , 

L-<L> = diag(ln(M2/Mi),0,ln(M2/M3)) . (58) 
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The mass M2 is considered as a kind of average of Mj. Usually, the 1st term of 
the right-hand side of Eq.(57) is considered. However, our model gives {YjY^)ij = 
(i 7^ j) because of the assumption (A. 3), and thus the LFV processes occur only 
through the 2nd term. 

In this model, the off-diagonal elements of m^jjCrriD become 

(m|)mD)i2 
(ml)mD)23 

Thus, T ^ and r ej processes are suppressed by factor {rncii/ 17102)'^ in 
comparison with the 67. While ^ ^ 67 is independent of Po, t ^ fi'j and 
T ^ 67 depend on (3q as well as cto- We find the branching ratios of £i £jj 
become larger as cos(q;o/2) become smaller. In the limit of cos(q;o/2) = 1, — > 67 
and T — > 67 do not occur, but t ^ /I'j can occur if /Sq is not or tt. 

For almost all values of cos(q!o/2) and /3o, we obtain 

Br(r /i7) < Br(r fxe) < Br(/i e^y) . (62) 

However, the branching ratios of these processes except for // — > 67 are too small 

to be observed in the future experiments. Their typical values are Br(£i — > ijj) < 
10-i2_ 



mm 



— m 



Dl 



a + b 
b 



In 



M2 



m 



Dl 



, Ml 
M3 



(59) 
(60) 
(61) 



4 Comments and discussions 

We constructed a model that all CP violation phases at the high energy scale and the 
low energy scale are controlled by the two Majorana phases which appear as the relative 
phases of neutrino masses. This strong restriction of the model is due to the assumption 
that Y^Yi, is diagonal. This requirement is motivated by the consideration that the large 
diagonal elements are needed to reconcile the maximal solar mixing angle at GUT scale 
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and the observed one at the low energy scale. We found that the induced l^isl is the 
measurable size in the near future experiments, the induced Dirac CP phase is likely 
in between — 7r/2 and — 37r/2, the asymmetry parameter can be of order 10^^, the LFV 
processes are suppressed. In this paper, we considered tan/3 = 20 case. For smaller 
tan/3, M3 becomes smaller as we see in Eq.(37). Then, Mi becomes smaller too. Since e 
is proportional to Mi in this model, e becomes smaller, so that it would become hard to 
explain the baryon number asymmetry in the universe. 

Of course, the assumption Y^Y^ may be too strong. The renormalization group argu- 
ment for the solar mixing parameters requires that the {YjY^)ii elements must be much 
larger than the other elements. Therefore, a general form would be 

/I * * \ 



yX = yl 



* * * 



\ * * * / 



(63) 



where y\ is of order 0.5 and elements shown by * are of order lO^"'. Although these 
elements are small, they will contribute to the LFV processes. However, their sizes are 
not controlled by the model and we lost the predictions for them. 

We can construct a similar model to the present model, by assuming the eigenvalues 
of the Dirac mass matrix are quasi-degenerate. This analysis will be reported soon. 
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Appendix A: Detailed derivations 

We parametrize as 

1 (Q 



where 5i are defined in Eq.(28) and 

Q 



mm \ 51^2 

1 
di 



(A.l) 



(A.2) 



We find by using a, b and c defined in Eq.(27), 

1 / QXQ SiS2QY 



ME' 



where 



a b\ I 



(A.3) 



(A.4) 



We block diagonahze the matrix M^^, by the seesaw calculation with respect to 
8182 <^ 1. We remind that a, 6, c are quantities of the same order because |mi| ~ \m2\ ~ 
Imsl- We find with 

1 -5^82Q-'X-^Y' 



8MQ~^x-^y)'' 



1 (QXQ 



4^1,1 \ 8f8lia + b-Y'^X-^Y) 



(A.5) 



(A.6) 



Next, we diagonahze QXQ in the first order of the small quantity 81 by applying 

/ 1 -8i{b/a) 0\ 



Vo 



8i{b/ay 

V 



1 
1/ 



(A.7) 



and we find that M^^ is diagonalized as 

{ViV2fME\ViV2)^dieig 
where Mj and phases are given in Eq.(31). 



' Qi<t>i ^i<f>2 g«</>3 ' 



Ml ' M2 ' M3 



(A.8) 
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